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Abstract 

The branch points of individual thermal self-energy diagrams at fc^ = 
4m'^,k'^ = 9m^, . . . are shown not to be branch points of the full thermal 
self-energy. Branch points of the full theory are determined by the complex, 
temperature-dependent energies of the quasiparticles, defined as the pole lo- 
cation, ko = £{k), of the exact retarded propagator. The full retarded self- 
energy is found to have branch points at kQ = 2£{k/2) and ko = 3£{k/3) as 
well as cuts in the space-like region. The discontinuities across the branch 
cuts are complex. The advanced self-energy is related by reflection to the 
retarded. 
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I. INTRODUCTION 



At finite temperature, self-energy functions have more branch cuts and more comphcated 
discontinuities than at zero temperature. The finite-temperature discontinuities have direct 
physical significance [1-9] . It is possible to compute the discontinuity of a self-energy diagram 
without having to compute the real part by employing cutting rules that replace certain 
propagators with Dirac delta functions [10-15]. All the known results about the location of 
branch cuts and the discontinuities across them apply at each order of perturbation theory. 
The pertubation theory is defined by choosing free thermal propagators that have poles at 
the zero-temperature mass m. This paper will demonstrate that when perturbation theory 
is summed the full self-energy will have branch cuts in different places and with different 
discontinuities than given in perturbation theory. 



A. Example at Zero- Temperature 

A simple zero-temperature example for a scalar field with interaction Cj = g4>^/6 will 
illustrate how higher order corrections can shift the location of branch cuts. Suppose that m 
is the physical mass but that one performs pertubative calculations using a free propagator 
A(A;) = 1/[A;^ — ml], where mo is some different mass. For simplicity mo should be finite 
and not the bare mass. The one-loop self-energy 

U^'\k)^'-^J0^^A{p)A{p-k) (1.1) 

has a branch cut for > 4ml. The discontinuity across the branch cut is 

Discn(i)(A;) = ^jd^p5+{p'-ml)5+{{p-kf-ml) 

[l-^f'\ (1.2) 
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The indication that k'^ — Ami is not a branch point of the full theory comes from the two- 
loop contribution. The full propagator is D'{k) = l/[k'^ — m^ — n(A:)] and by definition 11 
contains the necessary counter term to vanish at the true mass k"^ — m^. To do pertubation 



theory with mass mo the full propagator is written D'{k) — — ml — n(A;)] where 

n(A;) = — mo + Il{k). Of course IT does not vanish at fc^ = m^ or at fc^ = ml and this 
is the source of the problem. A self-energy insertion on the internal lines of (1.1) gives the 
two-loop self-energy term 

^^'\k)-ig'j ^[A(p)n(p)A(p)]A(p-fc). (1.3) 

This has a two-particle and a three-particle discontinuity. The quantity in square brack- 
ets has a double pole at — m^. Using [A(p)]^ = —d/S.{p)/dp^ lead to a two-particle 
discontinuity 

Discn(2)(A;) = ^jd^p 5;(p2 _ n(p) 

x5^{{p-kf-ml). (1.4) 
The presence of 5'(p^ — ml) requires an expansion of the shifted self-energy for p^ ^ ml: 

n(p^) = 5m^ + n(mg) + {/ - ml)Tl'{ml) + . . . 
where dm^ — m^ — ml- The integrated two-particle discontinuity is 



The second term changes the coefficient of (1.2) as required by wave function renormaliza- 
tion. The first terms is more important: It is infinite at — Ami- The infinity is a signal 

that the correct branch point is not at k"^ = 4ml. Multiple self-energy insertions on the 
same skeleton have two effects. First they modify the coefficient of (1 — 4mo/A:^)~^/^ to be 

5m' + U{ml) + Sm^ U'{ml) + -{5mYU"{ml) + ... 
^ 5m^ + U{m^) = 5m', 

where, in the last step, n(m') = has been used. Thus only Sm"^ survives as the coefficient of 
inverse square root. Multiple self-energy insertions also produce successively higher powers 
of the inverse square root: 
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1 2(5m^ 2.. 4mg 3/2 i 

This is the beginning of a Taylor series. All the corrections diverge at the false threshold 
k'^ = Ami- In the range Ami < k'^ < A{ml + 6m'^) each correction is finite but the Taylor 
series diverges. Thus pertubation theory fails throughout this region of k"^. To obtain a 
convergent series it is necessary to work in the range k^ > A{ml + Sm^). In this region the 
Taylor series converges and the sum is the full two-particle discontinuity with branch point 
shifted to the physical mass = ml + dm^: 

Discn(fc) = ^(l-^)V^ (1.5) 

The true two-particle threshold is still a square root branch point at k"^ = 4m^. The 
breakdown of perturbation theory is entirely due to a propagator A.{k) — l/[k'^ — ml] with 
the the wrong mass mo. The breakdown is easily avoided by using l/[k'^ — m^] for the free 
particle propagator. 



B. Non-zero Temperature 

In the previous example, individual diagrams of the perturbation series have branch 
points at the wrong threshold k"^ = Ami although the full theory does not. In finite- 
temperature field theory it is customary to perform perturbative calculations using free 
thermal propagators that have poles at the zero-temperature, physical mass m. With this 
choice the one-loop self-energy has a branch point at k"^ = Am^. However this is not a 
true branch point of the full theory. The insertion of the thermal self-energy on an internal 
propagator produces a two-loop correction analogous to (1.3) in which there is a double 
pole at = m^ because the one-loop self-energy does not vanish there. The double pole 
produces a discontinuity proportional to (1 — m^ /k'^)~^/'^ , which diverges at k"^ — Am^. This 
claim is easily checked by applying the Kobes-Semenoff cutting rules [10] to compute the 
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discontinuity. Both Le Bellac [14] and Gelis [15] display the two-loop discontinuity as an 
integral over c/^p containing 5'{p^ —w?). This is the same structure as in (1.4). It is computed 
explicity in Appendix A and the result is proportional to (1 — 4m^/A;^)~^/^, which is infinite 
at the false threshold just as in the T=0 example. 

The branch points of the full theory are not obtained by trivially replacing m? by a 
temperature-dependent effective mass. A proper calculation requires using unperturbed 
propagators with poles at the same energy at which the thermal self-energy vanishes so 
that there will be no double poles on internal lines. An energy £ which is a pole of the 
unperturbed propagator and also a zero of the self-energy is automatically a pole in the full 

— * 

propagator. Poles in the full propagator will occur at energy — 8{k) where E is is an 

— * 

complicated function of that depends on mass, coupling, and temperature. Moreover E 
is complex with the imaginary part being the damping rate of the single particle excitation. 
For definiteness the real part of E will be chosen positive and the imaginary part, negative. 
Thus E is in the fourth quadrant of the complex energy plane. The pole at — E{k) is 
called the quasiparticle pole. This paper will show that there is no branch cut at k"^ = Aw?. 
Instead there is a two-quasiparticle branch point in the full self-energy at the complex, 

— * 

temperature-dependent energy k^ — 2E{k/2). The branch point is the end point of a branch 
cut in which the two quasiparticles share the energy: /cq = E{ki) + E{k2) where k = ki + k2- 
The location of the quasiparticle pole in the full propagator is determined by effects 
that are higher order in the coupling. Approximating the full propagator by a simpler form 
that has a pole at the correct position reorders the perturbation series. This is similar to 
the Braaten-Pisarski resummmation [16,17] of high temperature gauge theories but differs 
in several respects. First, the breakdown of perturbation theory near the false thresholds 
is not an infrared effect. The breakdown occurs even in theories with masses and even if 
the temperature is small. Second, it is not necessary to retain the ko dependence of the 

— * 

self-energy in the new propagators, only the pole position E{k). 

A systematic method to organize the reordering of perturbation theory is to employ the 
integral equation that relates the full self-energy to the exact Minkowski propagator 
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and vertex F: 

Although one doesn't know the full propagator D^^, the natural first approximation is to 
use a free quasiparticle propagator that has a pole at the correct position. However the 
Minkowski integral equation is awkward to work with since it involves propagators with 
2^ components and vertices with 2^ components. It is simpler to use the imaginary-time 
formalism because there is only one propagator and one vertex function. In the imaginary- 
time approach the full self-energy is related to the full propagator T>' and vertex F by the 
single integral equation 

n(T, k)^^j dni2 j^dT'dr" V'{t', ki)V'{T", 4) 

xr{T',ki;T",k2;r,k), (1.6) 

where dQu — d^kid^k2S^{k — ki — fc2)/(27r)^. The existence of a quasiparticle pole at 
kg = £{k) in the Minkowski propagator determines the approximation to be used for the 
Euclidean propagator. By Fourier transforming n(T, k) and then analytically continuing, it is 
possible to obtain both the retarded and advanced self-energies TlR/A{kQ., k). This determines 
everything since each of the four real-time propagators D'^j^ik) are linear combinations of 
the retarded and advanced propagators. 

The paper is organized as follows. Sec 2 discusses the exact retarded thermal propagator 
D'^{k) and separates the quasiparticle pole from the self-energy effects. Sec 3 introduces the 
quasiparticle approximation to the propagator in both Minkowski and Euclidean space-time. 
Sec 4 computes the one-loop self-energy with quasiparticle propagators. The results are 
displayed in (4.5) to (4.8). The calculations are performed in the imaginary time formalism 
and then analytically continued to obtain 11^^/^. As a check, Appendix D performs the 
one-loop calculation entirely in the real-time formalism. The calculation is more difficult 
but produces exactly the same results. Sec 5 computes the two self-energy diagrams that 
contribute at two-loop order. The diagram in which there is a first-order self-energy insertion 
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on an internal line is the direct analogue of (1.3). Because the quasiparticle self-energy 
vanishes at ko = £, this diagram does not shift the location of the two quasiparticle branch 
point. The effect of this contribution is a only a change in the coefficient of the two- 
quasiparticle cut. Both two- loop diagrams have branch cuts for three-quasiparticle processes 
and these are computed. Sec 6 contains the conclusions and the general relation between 
the real-time Dab and Uab and the retarded/advanced quantities. 

II. EXACT PROPAGATORS 
A. Minkowski Space 

The quasiparticle poles occur in the Minkowski-space propagator and it is necessary to 
begin there and then convert to Euclidean propagators. The exact propagator in Minkowski- 
space has a 2 X 2 matrix structure. All four components are linear combinations of the 
exact retarded and advanced propagators D'j^{k) and D'^{k) as displayed in (6.2). Since 
^a{^) = ^r{~^) it suffices to investigate D'^{k). The retarded propagator is analytic in 
the upper-half of the complex /cq plane and satisfies the condition 

D'j,{k^,k)^[D'j,{-klk)]*- (2.1) 

At zero temperature the exact propagator has poles at ko — ±(m^ -|- k'^Y^'^. At non-zero 
temperature the location of these poles is temperature-dependent and complex. For defi- 
niteness, let the pole in the exact retarded propagator that occurs in the fourth quadrant 
be at ko — £ where 

E{k)^E{k)-iT{k)/2 (£;>0;r>0). (2.2) 

Both E and F are complicated functions of momentum, temperature, and coupling. The 
complex energy £ will be called the quasiparticle energy. Because of (2.1) the retarded 
propagator must also have a pole in the third quadrant at ko = —£*. Also because of (2.1) 
the residues of these two poles are related: 
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lim (A;o - S)D'ji{k) = Z/2E 

ko^S 

lim {ko + nD'nik) = -Z*/2E. 



(2.3) 



Here Z plays the role of the wave- function renormalization constant. The retarded propa- 
gator is directly related to the retarded self-energy 



D'j^{k) = k'^-m'^- Unik) 



(2.4) 



Because the full retarded propagator does not have poles at k"^ = rri^, the proper self-energy 
H-Rik) does not vanish at A;^ = w?. Therefore the usual Dyson-Schwinger expansion 



1 ^ ^R{k) 



+ 



nlik) 



+ 



(2.5) 



is not useful. The first term has a simple pole at k'^ = w?, the second term has a double 
pole, the third term has a triple pole,. . . . Performing perturbation theory around k"^ = w? 
is quite misleading. It is much better to write the full retarded propagator as 



Dnik) = \{k^ - S){ko + S*) - UR,p{k) 



-1 



(2.6) 



where the retarded quasiparticle self-energy is defined by 

Uji^pik) = Unik) +k^ + m^- \£\^ + tTko. 
By construction Ungp^k) vanishes at ko — S and also at ko = —£* 

Uj^p{S) = Tli^p(-S*) = 0. 
The natural expansion around the quasiparticle poles is 

1 'nRgp{k) 



(2.7) 



(2.8) 



+ 



+ ... 



(2.9) 



The second term has only simple poles at kg — S and at ko — —S*. It is convenient to define 

the derivative of the self-energy at these positions in terms of a complex constant B: 

(2EB ko^S 

(2.10) 

ko^-S* 



dko 



— < 



This constant B is related to the wave-function renormahzation constant in (2.3) by 



Z^l + B + B^ + ...^1/{1-B). 



(2.11) 



It will be helpful to have similar results for the advanced propagator. From the definition 



(2.12) 



the advanced propagator is analytic in the lower-half of the complex ko plane. It must have 
poles in upper-half plane at ko — £* and ko — —8. To emphasize these poles it is convenient 
to write the advanced propagator as 



D'^{k) = [(fco + S){ko - S*) - IlA,j>{k) 



(2.13) 



where the advanced quasiparticle self-energy is defined to be 



TlAqp{k) = TlA{k) + k'^ + m'^- |£:|^ - irA;o 



(2.14) 



B. Euclidean Space 

The finite-temperature Euclidean propagator is defined at discrete, imaginary frequencies 

where n is any integer. The full Euclidean propagator is 

( -D'(iLJn,k) ifn>0 
V'{iu;n,k) = { \^ - , (2.15) 

[ -D'j^{ioUn, k) if n < 
with the overall minus sign chosen for later convenience. Relation (2.12) for ko imaginary 
implies that D^(i27r|n|T, k) = D^(— i27r|n|T, k). It follows that V'^icUn, k) is an even function 
of n. The Euclidean propagator may be expressed in terms of the self-energy as 

T>'{iun, k) = {unf P Tl{iijjn, k) 

where the Euclidean self-energy is 



n(ia;„, k) 



IlR{iu!n, k) if n > 



(2.16) 



IlAii^n, k) if n < 

The relations Ilji{ko, k) = \U.R{-kl, k)]* and n^(A;o, k) = [n^(-A;o, k)]* guarantee that (2.16) 
is real. To emphasize the quasiparticle aspect the propagator may be written 



V'{iuJn, k) ^ - {i\uJn\ - S){i\uJn\ + £*) + Tlqp{iuJn, k) 



(2.17) 



where the quasiparticle self-energy is 



^qp{i^ni k) = Ii{iuJni k) + k"^ + — — T\u)r, 



(2.18) 



The presence of |n| rather that n in these results is very important but will cause complica- 
tions later. 



III. QUASIPARTICLE PROPAGATOR 

The natural approximation to the full Minkowski-space propagators is to retain the 
quasiparticle poles. Thus approximate (2.6) and (2.13) by 

Dnik) ^ 



{ko-S){ko + S*) 

The corresponding Euchdean propagator for free quasiparticles follows from (2.17): 

-1 



T>{iu;n, k) 



(z|cj„| - £){i\uJn\ + £*) 
1 

|2' 



Ujl + V\uJn\ + \S\ 

The transform to Euclidean time requires the Fourier summation 



(3.2) 



V{T,k)=T J2 e-''^"^V{iUn,k) (3.3) 

n=—oo 

for — /3 < T < p. Since (3.2) is an even function of the integer n, (3.3) is automatically an 
even function of r. To perform the summation it is convenient to write (3.2) without the 
absolute value bars on n as 
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Using this gives for the Fourier sum 



V{r,k) = A.([l + n(£)]e-^M + n(r)e^*H) 



2E 

This is the form of the quasiparticle propagator that will be used in the subsequent self- 
energy calculations. All the r dependence is of the form exp(— A|t|) where A is a member 
of the set below 

A e {S, —£*, —iuJi, —iw2, —i^z-i ■ ■ •} ImA < 0. (3.5) 
Each A has a negative imaginary part. The propagator will be written compactly as 

P(r)=E/(A)e-^'"' (3.6) 

A 

in which the coefficient functions are 

f\£) = [l+n{£)]/2E 
f{-S*) ^n{S*)/2E (3.7) 
f{-iue) = -2rra;,/(u;2 ^ g2^^^2 ^ ^*2y 

Although (3.4) will be used throughout, the infinite sum conceals several properties that 
are important to note. First, the time dependence exp(— A|t|) with Im A < will lead to a 
Euchdean self-energy that can be easily extended to the retarded self-energy in Minkowski 
space. However the starting point Vliuin, k) in (3.2) favors neither the retarded nor the 
advanced forms. Although it is not apparent, (3.4) is actually real: 

V{T,k)* = V{T,k). (3.8) 

This allows the time dependence to also be written exp(— A*|r|) if continuation to the 
advanced form of the Minkowski self-energy is desired. Although (3.8) is not obvious, it 
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must be true since T>{iijjn, k) is real and an even function of n. Appendix B proves (3.8) 
explicitly. Second, since exp(ia;„/3) = 1 the quasiparticle propagator (3.3) satisfies the KMS 
condition 

V{(3 -T,k) = V{T,k). (3.9) 

Without the infinite sum in (3.4) the KMS property would not hold. Appendix B proves (3.9) 
explicitly. Third, another way to obtain (3.4) is to begin with the time-ordered propagator 
in Minkowski space, which is given by the following linear combination of the retarded and 

advanced propagators: 

'''' - {ko -s){ko + s*) {ko + s){ko-s*y 

The Fourier transform, Dii{t, k), for real positive time t is determined by all the poles in the 
lower-half of the complex ko plane. These poles are at ko = S,ko = —S* and at ko = —icUn 
for n > 0. The propagator in Euclidean time results from continuing from positive, real t 

— * — * 

to negative, imaginary time —ir. The Euclidean propagator is 'D{t, k) = iL>ii(— ir, k) and 
gives precisely (3.4). 



IV. ONE-LOOP SELF-ENERGY 

It is always easy to perform loop corrections by integrating over Euclidean time and then 
Fourier transforming [1,18]. That method will be employed here. The first approximation to 
the integral equation (1.6) for the full self-energy is to use the quasiparticle propagator (3.4) 
and the bare vertex without corrections. This approximation treats the energy £ exactly 
even though it is a function of the coupling g. The one-loop correction shown in Fig. 1 is 

2 

n'(r, dnu V{t, hMr, h). (4.1) 

This may be expressed concisely using the notation (3.6) for the propagators: 

n'(T, k) = ^ [ dQu E /(Ai)/(A2)e-(^^+^^)l^l. (4.2) 

^ Ai,A2 
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The transform from r to discrete frequency a;„ is 

r/3 



U^iun, k)^ ( dre'^'-mHT, k) (4.3) 
= y / ^^^12 E /(Ai)/(A2)- 



n^^(A;o,A:) = I ycil]i2E/(Ai)/(A2) ^^_^^_^^ - (4.4) 



ia;„ - Ai - A2 

This can be extended from iujn for n > to complex ko with Im /cq > 0. It is analytic for 
Im /co > because Ai and A2 have negative imaginary parts. The extension therefore gives 
the retarded self-energy: 

,2 ^ 1 -(Ai+A2)/3 

/(Ai)/(A2) ^ 

A 

Although this is analytic for fco in the upper-half of the complex plane, when k^ is continued 
into the lower half-plane, the singularities at ko = h.\ + A2 produce branch cuts in the 
self-energy. 

Physical Cuts: It is useful to write out the various cases for the different Aj. First, if 
Ai = £1 and A2 = £2 the contribution to the self-energy is 

^2 d^i2 [l + n{£,)][i + n{£2)]-n{£,)n{£2) 



~2 J 2E{. 



(4.5) 



y2E2 ko-£i- £2 

The discontinuity across the cut is complex. The statistical factors provide for the Bose- 
Einstein enhanced emission of two quasiparticles minus the absorption of two quasiparticles. 
The second contribution is for Ai = £1 and A2 = —£2' 

r dQ,2 [1 + n{£,)]n{£;) - n{£,) [1 + n{£^)] 
2 J 2Ei2E2 ko-£i + £i ' ^ ' ^ 

The statistical factors account for a direct process in which quasiparticle 1 is emitted and 
quasiparticle 2 is absorbed minus the inverse process. If Ai = —£l and A2 = £2 the result is 



g- r d\h2 n{£l)[l + n{£2)\-[l + n{£l)\n{£2) 
2 J 2Ei2E2 ko + £t-£2 ■ 

If Ai = —£^ and A2 = —£2 the self-energy is 

g' r dQi2 n(£*)n(£*) - [1 + n{£*m + n{£^)] 
2 J 2Ei2E2 ko + £t + S^ ' 
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The one-loop self-energy is the sum of (4.5)-(4.8) and (4.9) displayed below. Appendix D 
computes the same quantity using the Minkowski propagators and obtains exactly the same 
answer. 

Unphysical Cuts: There are some additional contributions to (4.3). If both Ai and A2 
are positive integer multiples of — i27rT then the numerator of (4.3) vanishes. However if 
only one of the Aj is a positive integer multiple of — i27rT the numerator does not vanish. 

— * — * 

Since (4.3) is symmetric under interchange of ki <-> k2 it is only necessary to consider the 
case Ai = —ico^, A2 = 82 or —£2 and double the result to obtain 

2 fdnu^ 2TTiuJe 
9 / -^2^ 



1 1 

ko + iuji — £2 kQ + iue + £2 



X [ = I = 1 (4 9) 



These terms have branch cuts in the lower half-plane at kg — —iu!i+£2 and at kg — — io;^— 
The cuts are unphysical in that they are not entirely due to quasiparticle thresholds. The 
coefficient of this cut is proportional to the damping rate Fi and is in this sense a higher 
order effect. Sec 5 will show that (4.9) is exactly canceled by two-loop effects. For later 
comparison it is useful to return to the term in H (r, k) whose frequency transform produced 
this cut: 



. 00 
g'J rfJlisEe'"^' 



2T Tiui 

X [[1 + ^(^2)]^e-^- + n{£;)^S^] . (4.10) 
Advanced Self-Energy: Since quasiparticle propagator satisfies the KMS condition, the 

— * — * 

integrand of (4.1) could equally be written D(/3 — r, ki)T>{j3 — r, k2). The Fourier transform 
to iuin is then expressed as 

-,2 , ^-(Ai+A2)/3 _ I 



u.\iujn,k) = I dni2 E /(Ai)(A2)- 



Ai,A2 + Ai + A2 

This is exactly the same self-energy as (4.3). However in this form it is easily extended from 
to a function of complex kg that is analytic for Im /cq < 0. This extension gives the 
advanced self-energy 
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n^(fco,fc) = ^ / rf^^l2 E /(Al)/(A2) , ^ , ^ , ■ 
^ Ai,A2 /Co + Al + A2 

It satisfies n^(A;) = ni<;(— /c) as required and has all its branch points in the upper half of 
the complex kg plane. 

Mixed Representations: Because of the KMS condition one can also represent the self- 

— * — * 

energy using a mixed form V{t, ki)V{P — r, k2). This leads to 
U\iu;n,k) = ^ / dn^, E /(Ai)(A2). 

I J ^^^^^ ILOn — Ai + A2 

Although this is the same self-energy, this representation cannot be easily extended to either 
the retarded or the advanced form of the self-energy. In Sec 5B it will be necessary to use 
the KMS identity in a similar way to manipulate the two-loop self-energy into a form whose 
Fourier transform will be analytic in the lower half-plane. 

V. TWO-LOOP SELF-ENERGY 

The simplicity of the one-loop calculation makes it likely that the two-loop contributions 
can be computed by the same method. The contributions of Figs. 2 and 3 will be denoted 
by and 11^^ respectively. 

A. Self-Energy Insertion on Quasiparticle Propagator 

The value of the diagram shown in Fig. 2 is 

Ha (r, k) = -g^ J d^u V\t, hWr, h), (5.1) 
where is the one-loop corrected propagator: 

V\t, k) = fdT'dT"V{T - r", A^)n,p(r" - r', ^)P(r', k). 

J 

This is not the most convenient way to compute T>^ . It is easier to employ the method 
discussed after (3.10). This requires the Minkowski-space time-ordered propagator, now 
with one insertion of the retarded and advanced self-energies: 
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{ko - sY{ko + s*Y {ko + s)^iko - s*Y • 

To Fourier transform to real, positive time t requires closing the ko contour in the lower-half 
of the complex k^ plane. The singularities in ko in the lower-half plane are as follows: (1) a 
simple pole at /cq — ^i (2) a simple pole at k^ — —£*, (3) simple poles in n(/co) at /cq — —i^i-i 
and (4) branch cuts in Iljiqp{k). Consequently the Fourier transform is 



iD[,{t,k) = [I + n{£)]^e-^'' 



R* 



{^{k^-8Y{k^ + £*y 



ko=—iuJi 



{^{ko + £y{,ko-£*y 



ko=—iu)i 



/-cute dk.j . 
J/m fcn<0 In 



■ikot 



The self-energies IlRqp and H^gp can be expressed in terms of Ur and 11^ using the definitions 
(2.7) and (2.14). Evaluating the propagator at the Euclidean time t — —ir gives 

V\r, k) = [l + n{£)]^e-'^ + n{r)^e''^ 



,tU>(T 



1=1 



m 



88*) + T^ujj] 



(^2 + £2)2(^2 + ^*2)2 

^^(-^a;^) n|j(ia;^) 



{ujj + 88* -Vuj^Y {ujj + 88* + Vujiy 



^ routs 

Jim ko<0 ^TT 



(5.2) 



One way of proceeding is to add this correction to the free quasiparticle propagator (3.4). In 
the sum T> + T>^ the coefficients of the quasiparticle terms are modified to 1 -|- S and 1 + B* 
and the term proportional to F in (3.4) cancels in the sum. It was this term that produced 
the unphysical cuts in the one-loop self-energy. The cancellation inV + guarantees that 
unphysical one-loop cuts will be canceled in two-loop order. The following discussion shows 
these features explicitly as well as the three-quasiparticle cuts that arise. 
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Wave Function Correction to the Two Quasiparticle Cut: To compute the self-energy 
requires substituting (5.2) into (5.1). The contribution of the first hne of (5.2) to n^^(T, k) 
is 

^1 „-£ir , „/c7*N ^1 



-g'J dn,,[[l + n(£i)]^e-^- + n{S*,)^e'^ 



X 



2Ei ' ''2Ei 

[1 + ^(^2)]^e-^- + ^(^2)^e^^] . (5.3) 



This may be symmetrized so that Bi and B2 appear equally. When added to (4.2) it merely 



introduces the wave function correction 1 — Bj ^ Zj. 



Cancellation of Unphysical Cuts: The r dependence in (4.10) produced the unphysical 
cuts in (4.9). When the second line of (5.2) is substituted into (5.1) it gives 



/oo 
dQi2^e 
1 



X [[1 + ^(^2)]^e-^- + n{El)^S^\ . (5.4) 

This exactly cancels (4.10) so that the one- loop unphysical cuts are removed. Obviously the 
third and fourth lines of (5.2) will produce new unphysical cuts in the two-loop self-energy. 
These will be canceled by higher loop effects. 

Cut for Three Quasiparticles: The last term in (5.2) requires integrating in ko around 
the branch cuts in the one-loop self energy: 

• dko n^(/c)[l + n(fco)] 

Jim fco<o 2n {ko - SiYiko + SlY ' ^ ' ^ 

It is convenient use the representation (4.4) but to change the internal momentum variables 
to ks and k^ in correspondence with Fig. 2: 

1 _ g-(A3+A4)/3 

/(A3j/(A4j 

A3,A4 



n^(^o) = C [dQs, Y: /(A3)/(A4)V- 



A3 - A4 



The denominator ko — A3 — A4 produces the branch cut in ko- The integration around the 
cut is performed by interchanging the order of integration to get 



^ f,^ ^ /(A3)/(A4)e-(^3+A4)r 

2 J '\ti, (A3 + A4 - s^nAs + A4 + Sir ■ 
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This is the exphcit evaluation of T>1^^{t, k), i.e. the last hne of (5.2). When substituted into 

^ — * 

(5.1) the contribution to 11^ (r, /c) is 



- /(A2)/(A3)/(A4)e-(^^+^3+A4)r 
(A3 + A4 - ^i)2(A3 + A4 + £t)^ ■ 



This is easily transformed to get n^^(ia;„, A;). The extension from icun to complex ko analytic 
in the upper half-plane is 

/(A2)/(A3)/(A4) 



{A} 

g -(A2+A3+A4)/3 _ ]_ 
(A3 + A4-£l)2(A3 + A4 + £1*) 



A4 



This contains the cuts for three quasiparticles at /cq = A2 + A3 + A4. The unphysical values 
of A will be canceled by higher loops. This completes the analysis of Fig. 2. 



B. Vertex Correction 

Fig. 3 shows the two-loop diagram containing a vertex correction. Two of the loop 

— * 

momenta are independent. For definiteness the independent momenta are taken as ki and 
^3 and dn = d^k^d^k^/ {2t:)^ . The remaining k2,k4,k5 are linear combinations of the these 

— * 

two and the external k. The self-energy contribution is 

^Bir,k) fdQ f^dr' dr" V,ir')V2ir") 

A J Jo Jo 

xV,{t" - t)V,{t' - t)V,{t" - r'). (5.7) 

The three times r, r', and r" lie in the interval [0,/3] and may be ordered in six different 
ways as follows: 

B1:t' < t" <t B2:t" <t' <t 
B3:t<t'<t" BA:t < t" < t' 
Bb:T'<T< t" B6:t" <t< t' 
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The left and right columns differ by an interchange of r' and t" . Because of the structure of 
the integral, this is the same as interchanging Ai ^ A2 and A3 <-> A4. Thus only Bl, B3, and 
B5 need to be computed. With the representation (3.6) for the quasiparticle propagators, 
the integration over Bl gives 

En 

{A}i=l 

X 



g-(Ai+A2)r 



(Ai - A4 - A5)(Ai + A2 - A3 - A4) 



+ 



g-(A3+A4)r 



(A2 - A3 + A5)(Ai + A2 - A3 - A4) 



+ 



g-(A2+A4+A5)r 



(5.8) 



(Ai - A4 - A5)(-A2 + A3 - A5)- 

The T dependence of these three terms will easily lead to two-particle cuts at /cq = Ai + A2, 
/co = A3 + A4, and a three-particle cut at /cq = A2 -|- A4 -|- A5. The next integration, B2, gives 
the same answer as (5.8) but with the interchanges Ai ^ A2 and A3 ^ A4. 

Integration B3 can best be done by using the KMS condition to rewrite it as 

^m{rA = jjdn J^t' £ dr" V,{t')V,{t") 
xVsiP + T- t")V^{(5 + t - t')V^{(5 + t' - t"). 

The time argument for each of the quasiparticle propagators is positive. For example, for 
the time dependence is exp[— A3(/3 + t — t")]. The integrand written in this form leads 
to the most convenient form for the final answer with n(T) a product of terms of the form 
exp(— At) as desired. Direct integration gives 

,4 ^ 5 



nM.(r,^) = ^/rff^En/(A.) 



X 



{A}i=l 

g-(Ai+A2)T g-(A3+A4+Ar,)a 

(A2 - A3 - A5)(Ai + A2 - A3 - A4) 

g-(A3+A4)r g-(Ai+A2+A5)/3 



(Ai - A4 + A5)(Ai + A2 - A3 - A, 

t 



g-(Ai+A3+A5)r g-(A2+A4)/3 



(Ai - A4 - A5)(-A2 + A3 + A5) 
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(5.9) 



The tau dependence of these terms will again produce two particle cuts at ko — Ai + A2, 
kg = A3 + A4, but a different three-particle cut at ko — Ai + A^ + A^. Integration B4 requires 
interchanging Ai A2 and A3 A4. 

The contribution of B5 is more difficult. First use the KMS condition to write it as 

xP3(/3 + T - t")V^{t - t')V^{I3 + t' - t"). 

The integration gives 

K,{r,k)=^^jdnY.\lf{A,)C 



X 



{A}j=l 

g-(Ai+A2)r g-(A3+A5)/9 _^ g-(A3+A4)T ^-Aa/? 



_g-(Ai+A3+A5)T g-A2/3 _ g(-A2-A4+A5)r ^-(As+As)/? 



^ - (Ai - A4 + A5)(-A2 + A3 + As) ■ ^^-^^^ 
The last term contains tau dependence exp(-|-A5T) which, when Fourier transformed, is 

difficult to extend analytically in the upper half-plane. It is useful to isolate all the A5 

dependence of this term by defining 

A5(/3-r) 

Q = . 

(Ai - A4 + A5) (- A2 + A3 + As) 

The generalized KMS relation (C9) proven in Appendix C shows that 



As As - (Ai - A4 - As) (-A2 + A3 - As) 

g(-Ai+A4)r^^(_^^ + A4) - e(-^2+A3)(/3-r)^^(_^2 ^ 



Ai + A2 - A3 - A4 

where F is the function defined in (CI). When this is substituted into (5.10) the result is 

{A}j=l 

g-(Ai+A2)r g-(A3+AB)/3 _^ g-(A3+A4)r g-A2/3 _ g-(Ai+A3+A6)T g-A2/3 

-^^^'^ (Ai-A4 + A5)(-A2 + A3 + A5) 

g-(A2+A4+A5)T Q-AsP 

•^^'^'^(Ai-A4-As)(A2-A3 + A5) 

^g-(Ai+A2)r e-^^^F^i-Ai + A4) - e-(^3+A4)r g-A2/3^g(_^2 ^ 

Ai A2 - A3 - A4 
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(5.11) 



The T dependence determines the ko dependence. The terms exp(— (Ai + A2)t) and 
exp(— (A3 + A.4)t) produce two particle cuts at ko — Ai + A2 and A;o = A3 + A4. The 
terms exp(— (Ai + A3 + A5)t) and exp(— (A2 + A4 + A5)r) produce three particle cuts at 
/co = Ai + A3 + A5 and /cq = A2 + A4 + A5. Integration B6 requires interchanging Ai ^ A2 
and A3 <-> A4. 

VI. CONCLUSION 

The above results follow from the existence of poles in the full retarded propagator 

D'j^{ko, kj) at energies ko = Xj where 

Xj = £(kj) or - S*{kj) ImA < 0. (6.1) 

These poles were shown to produce singularities in retarded self-energy integrands. In the 
two-quasiparticle channels there are singularities at /cq = Ai + A2. In the three-quasiparticle 
channels the singularities are at kg — Ai + A2 + A3. Contributions with X — S correspond 
to stimulated emission of quasiparticles weighted by 1 + n{£); contributions with A = —£* 
correspond to absorption of quasiparticles weighted by n{£). 

The singularities in the integrands of IlR{k) produce branch points when they are trapped 
at end points of the three-momentum integrations. Without knowing the momentum de- 
pendence of £{k) it is only possible to analyze this trapping in the equal mass case, i.e. 
when all the internal lines have the same dispersion relation £(k). In that case the pole of 
the integrand at ko — £{ki) + £{k2) produces an end point singularity from ki — k2 — k/2. 
The branch point is thus at ko = 2£{k/2). For 3 quasiparticles the branch point is at 
ko = 3£^(A;/3). The poles of the integrand at ko = £{ki) - £{k2)* and kg = —£{ki)* + £{k2) 
produce end point singularities from the region k\ — ak, k2 — {1 — a)k where a ±00. 
Since all radiative corrections vanish at infinite momentum, the branch points are near the 
real axis at ko — ±\k\ — irj. These results hold only for equal masses. In general the branch 
point locations will depend upon the functions £{k). 
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Cuts in the retarded propagator automatically give those of the advanced propagator 
because D'j^{k) = D'j^{—k). This also determines the four real-time propagators 

D[,{k) = [1 + n{ko)]D'^{k) - n{ko)D'^{k) 
D',,{k)=e^^-'^^'^n{ko)[D'^{k)-D'Ak)] 

D'^^ik) = n{ko)D'j,{k) - [1 + n{ko)]D'^{k). (6.2) 

Each branch cut of the is completely below the real axis or completely above. There 
are no branch cuts that cross the real axis. In addition the have simple poles at 
ko = ±i27rnT from the Bose-Einstein functions. Although the D'^f^ can be written in terms 
of the the thermal Feynman propagators Dp^p this introduces step functions 9{ko) which 
make the analytic properties of Dp^p more complicated. 

The real-time self-energies are related to the inverse full propagator by 

[D'{k)]-J = ie-m')al,-U^,{k). (6.3) 

In terms of the retarded and advanced self-energies this imphes 

Unik) = [1 + niko)]Un{k) - n(A;o)n^(A;) 
ni2{k) = e'''''>n{ko)[-nR{k) + UA{k)] 
U^iik) = e(^-'^)'=°n(A;o)[-n^(A;) + n^(fc)] 

U22{k) = n{ko)UR{k) - [1 + n{ko)]UA{k). (6.4) 

Several interesting points require further investigation. The separation of free quasipar- 
ticle effects was done by rearranging the propagator. It would be useful to have a operator 
method for separating the free quasiparticles from the interactions. Work on this is in 
progress. A related problem is whether the discontinuities can be computed directly with- 
out having to compute the entire self-energy as done here. In the perturbative approach, 
the cutting rules of Kobes and Semenoff [10] accomplish this. However their derivation also 
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requires using the operator structure. The physical significance of the discontinuities re- 
quires further investigation. Since the true branch points he off the real ko axis it is natural 
that the discontinuities across the branch cuts are complex. For example, the two-particle 
discontinuity of (4.5) is 

Discn^(^) = -z^ / ^f^27r5(^o - ^1 - S,) 

x[[l + n{£,)][l + n{£2)] n{£,)n{£2)]. (6.5) 

This is very much like what would be expected for the difference between the production 
rate of two quasiparticles minus their absorption rate, except that the quasiparticle energies 
S are complex. 
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APPENDIX A: BREAKDOWN OF PERTURBATION THEORY 

If one applies the Kobes-Semenoff cutting rules [10] to Fig. 2 using free thermal propa- 
gators it has the same breakdown near threshold as the T — example discussed in Sec lA. 
The formula for this particular discontinuity is displayed in Le Bellac [14] and in Gelis [15]. 
The two-particle discontinuity is 

Discni^(A;) = J d'^p [1 + n(po) + n{ko - po)] 
xe(po)(^'(p'-m2)RenR(po) 

xe{ko-po)S{{p-kf-m^). (Al) 

The contribution of Imll/j has been dropped since it produces a three-particle discontinuity. 
To display the result it is useful to let k — \k\ and — k^ — and a — {1 — 4m^/i^^)^/^. 
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Direct integration gives 



-ig^ kp + ak kp - ak 



X 



Jco — a/c 
(A; + Q;fco)RenR( 



+(k - ako)Ren.R{^^ ) e(ko - ak) (A2) 

where kinematics requires that either K"^ < or > 4m^. At the perturbative two-particle 
threshold, K'^ — > 4m^, so that a — > and 

^^'""^^^^ ^ 16^[^ + 2n(|)]Ren^(^) (A3) 
The behavior of this discontinuity like (1 — Am^ j K'^)~^l'^ produces an infinite correction 
at the false threshold which signals the breakdown of perturbation theory just as in the 
zero-temperature example of Sec lA. One can also check from (A2) that at the lightcone 
threshold, K"^ —>■ 0~, the discontinuity does not diverge. In retrospect, this is because the 
quasiparticle effects do not change the location of the space-like branch cut for equal masses, 

— * — * 

— < /co < 1^; as discussed in Sec 6. 

APPENDIX B: REALITY AND KMS CONDITIONS 

It is not obvious that the quasiparticle propagator 'D{t, k) displayed in (3.4) and used 
throughout the paper satisfies the reality and KMS conditions claimed in (3.8) and (3.9). 
The infinite sum in (3.4) obscures these properties. One can rewrite that sum in another 
way using 

gia;„|r| ^ g-ia;„|r| ^ 2i sin(a;„|T|). 

The sum over sin(a;„|T|) can be performed using the identity 

im{uJn\r\)j^^ 

n=l 



T^sin(u;„|r|)— 2 , c*2^ 



^ -[l + n(£:)]e-^l^l -n(r)e^*l^l 



4:E 



+[1 + n(£*)]e-^*l^l + n(£)e^l^l] . (Bl) 
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Using this in (3.4) gives 

P(r, k)^^ [[1 + n(£:*)]e-^>l + n{S)e'^^\ 

Each term on the right hand side is the complex conjugate of the corresponding term in the 
original expression (3.4). Hence V{T,k) is real. 

To prove that the quasiparticle propagator satisfies the KMS condition requires writing 
the propagator in yet another way. In the original form (3.4) use 

gia;n|T| ^ cos(a;„|T|) + i sin(a;„|T|). 

The sum over sin(a;„|r|) can be performed with the identity (Bl) to give the result 

D(r, k) = ^ \[l + n(£)]e-^l^l + n{8*)e^*\^^ 
+[l + n(£:*)]e-^*l^l+n(£:)e^l^l" 

oo 'iVoj 

In this form the KMS condition V{(3 — t, k) = V^t, k) is satisfied manifestly. 



APPENDIX C: GENERALIZED KMS IDENTITIES 

In Sec 5B it is necessary to use some relations that are generalizations of the KMS 
identity. To demonstrate these it is useful to define 

^ {k^-£){k^ + S*) {ko-£){ko + S*y ^ ' 

This satisfies 

F(-fco) = e^'=°/(^o)- (C2) 

F has poles in the lower-half of the complex ko plane at ko = A where A e {S, —S*, —iujn}- 
At the poles 
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where /(A) are the functions given in (3.7). It also has poles in the upper half-plane at 
ko = -A: 

KMS Identity: Because F{ko) vanishes sufficiently rapidly in all directions of the complex 
plane as \ko\ oo, the contour integral (C3) vanishes when the contour C is a circle of 
infinite radius: 

0=1 ^ F{ko) e^'^^-^^ (0 < T < 13). (C3) 

The vanishing of the integral imphes that the residues of the lower half-plane poles cancel 
those of the upper half-plane: 

^/(A)e-^(^--) = ^/(A)e-^^ (C4) 

A A 

Since the left and right sides of this are the Euclidean propagator (3.6), this just proves the 
KMS theorem 

V{f3 -T,k) =V{T,k). (C5) 

Theorem 1: For C a circular contour at infinity and x any complex number inside the 
contour, the following integral vanishes 

r dkn pkoW-T) 

The contribution to the integral of the poles at ko — A, ko — —A, and ko — x must all cancel. 
This implies 

This is a generalization of the KMS identity. If the differential operator [x + d/dr) is applied 
to both sides of (C7) it reduces to (C4). 
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Theorem 2: For the same contour as before and < r < /3 the integral (C8) vanishes 
Evaluating the integral by Cauchy's theorem gives 

g-A(/3-r) g-Ar 

5 ^^^^ (A + x)(A + |/) ^ ? ^^"^^ (-A + a;)(-A + 2/) 

Applying {y + d/dr) th this reproduces (C7). This identity is used in rewriting (5.10) in the 
form (5.11). Obviously these identities could be genralized to polynomial denominators of 
any order. 



APPENDIX D: ONE-LOOP CALCULATION IN THE REAL-TIME FORMALISM 

Calculations may also be done directly in the real-time formalism. This Appendix will 
compute the one-loop self-energy in the real-time formalism and show that the answer is 
the same as obtained rather easily in Sec 4. In the quasiparticle approximation the real- 
time propagators Dah{k) are the linear combinations (6.2) of the approximate retarded and 
advanced quasiparticle propagators 



{ko-S){k^ + S*) ' {k^ + S){k^-S*) 

The retarded self-energy that implied by (6.4) is 

{e^^° + 1) Unik) = e^^" nn(A;) - Tl22{k). (Dl) 

The one-loop contribution has two propagators with momenta ki and ^2 . Integration will 
be over ki with the other defined hy k2 = ki — k. The necessary one-loop self-energies are 

ig"^ f d^ki 

" " " v2) 



nn(^) = \j j^^n{k,)n{k2)[e^''WR{k,) - D A{k,)][eP'^W nik^) - D^ik^ 
ll22{k) = ^ I ^^n{k,)n{k2)[Dn{k,) - e^^^^ D A{k,)][D nik^) - ef^'-WAik^)]. (D2) 
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When these are substituted into (Dl) the term DA{ki)DR{k2) cancels. The remaining three 
products of the form D{ki)D{k2) are multiphed by combinations of exponentials that cancel 
one of the Bose- Einstein functions n{ki) or n{k2). The result is 

(e^^° + 1) Un{k) = ^ I ^^n{k,)DA{k2W'' + l)Dn{k^) - {e^''° + e'^^-)D^(^i)] 

+1" / - (e^^° + l)DA{k2) + {e^'^^ + l)Dn{k,)]. (D3) 

Note that D^{ki)DA(k2) appears in both lines. It is convenient to compute the first integral 
by closing the /cqi contour below. The poles in the lower half of the /cqi come from two 
sources: Dji{ki) has quasiparticle poles at koi — £i and koi — —£*] and n{ki) has poles at 
/coi = —iuj^. After the fcoi integration is performed, there is a common factor e'^*^" + 1 on the 
right hand side. The contribution to from the first line of (D3) is 

n{-£l) 



f d^ki 1 n{£i) 

2 J (27r)3 2Ei ^ (ko -Ex- £2) {ko - Si + £2) {ko + El - £2) (ko + £{ ^ ElY 

V ffk^r^Y^ I 2iu;e{£l-£,) 

^2 J {2nY ^ (A;o + iue - £2) {ko + icoe + <?|) i^e + ^1) i^e + ^f) ' 



(D4) 



where 2Ei = £1 + £1. To compute the integral on the second line of (D3) it is convenient to 
close the koi contour above. The poles in the upper-half of the koi plane come from D^(A;2) 
(recall ko2 — koi — ko) at koi — /cq + ^I ^01 — ko — £2 and from the Bose- Einstein function 
n{k2) at koi = ko + iuie- The second line of (D3) contributes 



d^ki 1 ^ n{£*) n{-£2) 



(27r)3 2E2 ' {ko + £i - £i){ko + £i + £*i) {ko -£2- £i){ko -£2 + £*i) ' 



+^/^TV I 2zo;,(g* - £2) 

^ 2 J (27r)3 (^0 + iooe - ^1) (^0 + iooe + Sf) {uj + £|) + £f) ' ^ > 

The sum of (D4) and (D5) gives Il.R{k) to one-loop order. It agrees completely with the sum 
of (4.5)-(4.9). In this method of calculating, the unphysical branch cuts produced by the 
denominators containing ko + ioJe. + z arise from poles in the Bose-Einstein functions. They 
are not artifacts of the Euclidean calculation performed in Sec 4. 
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Fig 1: One- loop self-energy. 



Fig 2: Two-loop self-energy due to one self-energy insertion. 



Fig 3: Two-loop self-energy due to vertex correction. 
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